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Abstract 

Many-particle simulations of vehicle interactions have been quite successful in the 
qualitative reproduction of observed traffic patterns. However, the assumed interac- 
tions could not be measured, as human interactions are hard to quantify compared 
to interactions in physical and chemical systems. We show that progress can be 
made by generalizing a method from equilibrium statistical physics we learned from 
random matrix theory. It allows one to determine the interaction potential via dis- 
tributions of the netto distances s of vehicles. Assuming power-law interactions, we 
find that driver behavior can be approximated by a forwardly directed 1/s poten- 
tial in congested traffic, while interactions in free traffic are characterized by an 
exponent of a ~ 4. This is relevant for traffic simulations and the assessment of 
telematic systems. 
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1 Introduction 



In recent years, statistical physics and nonlinear dynamics have been very 
successful in modeling, simulating, and understanding empirically observed 
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instability and pattern formation phenomena in freeway traffic [1,2,3,4,5,6]. 
On the one hand, vehicle traffic flow can be well approximated as a particular 
one-dimensional particle gas with Boltzmann statistics [6,2]. On the other 
hand, the motion of a vehicle i of mass m and length l { at location rj(i) with 
speed Vi(t) = dri/dt, maximum velocity t>o, adaptation time r, and fluctuation 
force £i(t) can be reflected by an acceleration equation of the form 

dvi v -Vi 

m—=m + f{si)+U t )- I 1 ) 

dt r 

f(si) is a negative and repulsive force, which depends on the netto dis- 
tance (clearance) = (r i+1 — r\ — /j) between two successive vehicles. With 
Wi(s) = v o + r/(s)/m, Eq. (1) corresponds to the optimal velocity model 
dvi/dt = [Wi(si) — Vi]/r [7] with an additional noise term £j. Due to a lack 
of good single-vehicle data and suitable evaluation methods, the form of the 
interaction potential U(s) with f(s) = —dU(s)/dr = —dU(s)/ds has been 
subject to speculation. A direct evaluation by averaging v j-over-Sj data suffers 
from the wide distribution of netto distances and causes problems of interpre- 
tation [2,8]. Exploiting the source of this problem, we suggest a statistical ap- 
proach to determine the potential and present first results regarding its form, 
based on a comparison of single vehicle data with the netto distance distribu- 
tions obtained from rigorous solutions for a particle gas in equilibrium. For a 
justification of this approach see Ref. [9]. 

Normally, the interaction forces or potentials between particles governing each 
other's motion, are not directly measurable. The great success of scattering 
theory was to determine interaction potentials from statistical distributions of 
particles scattered at some "target" composed of the material under investi- 
gation [10]. Learning about human interactions requires a somewhat different 
approach from statistical physics, which can be learned from books on ran- 
dom matrix theory [11]. This method was successfully applied to the statistical 
description of the time gaps Tj between the arrival times of buses in some Mex- 
ican cities [12], where the potential U was found to be a logarithmic function, 
i.e. U(Ti) = — ln(Tj). It was speculated that the same potential would approx- 
imate the (time) headway distribution of highway traffic [13]. In contrast, our 
study will identify the shape of the spatial interaction potential and reveal a 
different driver behavior in free and congested traffic. 



2 Methodological Approach 

We have recently extended a method developed for classical many-particle 
systems exposed to a "thermal bath" of a given temperature, i.e. to random 
forces of a certain variance and statistics [9] . The resulting velocity and netto 
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distance distributions [see Eqs. (7) and (10)] allow one to draw conclusions 
about the interaction potential U, as this determines their shapes. These dis- 
tributions have been originally derived assuming a conservation of momentum 
and energy, i.e. a transformation of potential energy into kinetic energy. How- 
ever, by investigating a Fokker-Planck equation equivalent to Eq. (1), one can 
show that the same distributions are steady-state solutions for driven many- 
particle systems with forwardly directed rather than symmetrical potentials 
[9], if the system is large enough and if the average velocity V and variance 9 
are constant. For dV/dt 0, d6/dt ps 0, and iV 3> 1, fluctuations in the sys- 
tem become negligible, and the total energy is only slightly fluctuating. The 
effects of the driving force mvo/r and the dissipative force —mvi/r balance 
each other in a statistical sense [9]. However, dV/dt ~ requires linear stabil- 
ity of the system (1) of coupled differential equations, i.e. dW(s)/ds < l/(2r) 
[7]. Therefore, we have to exclude the stop-and-go regime between 20 and 40 
vehicles per kilometer and lane from our investigation [2]. In summary, the 
data analysis must be carried out for ensembles of vehicles in a stationary 
traffic state with a well-defined velocity variance (generalized "temperature" ) 
Var(-Uj) = {crV) 2 in a reference frame moving with a constant velocity, the 
average vehicle velocity V. That is, one has to restrict to small density inter- 
vals, otherwise one will mix up systems of different generalized temperatures 
and different average velocities. Such a careful analysis is presented here for 
the first time for single vehicle data of the Dutch two-lane freeway A9, which 
was confirmed by another analysis for Czech highway data (not shown). This 
has led to new and surprising insights regarding the shape of the interaction 
potential and its dependence on the traffic state. We will focus on forwardly 
directed power-law potentials 



U(s) oc < 



s a for s > 
otherwise , 



(2) 



where a > is a fit parameter and s the netto distance (clearance) between 
two successive cars on a (ring) road of length L. Similar relations have been 
suggested in some car-following models [14,15] and describe the repulsive ten- 
dency of drivers to keep a safe distance to the respective car ahead. The 
exponent determines the characteristic dynamic and stationary behavior. 



3 



3 Short-ranged Dyson's gas with power-law potential 



We will now investigate the statistical gas of iV point-like identical particles 
on a ring of scaled length N interacting via the potential energy [16] 

N 

U(xi,...,x N ) = C^2U(x i+1 -Xi). (3) 
i=i 

Herein, (xi, . . . ,xn) is the vector of scaled particle positions 

i-l jy i-l 
3=1 L 3=1 

where p denotes the global vehicle density. For convenience, we will use a pe- 
riodic index, i.e. x i+ ^ = Xi + N. Let this so-called short-ranged power- law 
Dyson's gas (SRDG) be exposed to a temperature reservoir with Boltzmann 
statistics and generalized temperature a 2 > 0. The corresponding dimension- 
less Hamiltonian reads 

i N C N 

with the additional condition Y^f=i \ x i+i ~ x i\ — N, where Ui = Vi(u) /V is the 
scaled velocity of the ith particle and (u) = 1 the scaled average velocity. Note 
that all the quantities used here are dimensionless. The probability of finding 
the system in the phase-space element Q = (xi, . . . , xn, u±, . . . ,un) is 



N 



i=i 



N / N \ 

x -q e -M*i + i-*i) s In - \ x i+i - *i\ , (6) 

i=l V i=l ) 



where (3 = C/{ma 2 V 2 ) represents the scaled generalized inverse temperature 
and M the normalization factor, which can be determined from the normal- 
ization condition J r2 n P* (fi) dVt = 1 [16]. By integration over the 2N — 1 
independent variables, the probability of finding the ith particle with velocity 
Ui is obtained as 



(^-{u)f 
2a 2 

That is, the velocities are Gaussian-distributed with mean value 



p \ u i) = -7==- exp 
\ / 2na 



(7) 



(u) = I uP'iu) du = 1 (8) 
Jr 



4 



and variance 

( {u - (u)f) = f(u- (u)) 2 P'(u) du = a 2 = ((v - V) 2 )/V 2 . (9) 

JR 

For the above SRDG, one can also determine the equilibrium distribution 
of scaled particle distances Zi = (x i+1 — Xj) = pSi [17], which gives us the 
expression for the netto distance (or clearance) distributions: 

P^\z)=Ae' fSz ' a e- Bz . (10) 

Herein, A = A(a,(3) and B = B(a,j3) are constants determined via the nor- 
malization conditions 

j™pf ) {z)dz = 1, (z)= J™ zPj?\z)dz = l. (11) 

These equations can be analytically solved only for particular potentials, in- 
cluding power-law relations. In the following, we need two special variants of 
the SRDG. In the case a — 1, we find 

J~P$\z)dz = 2A^JC 1 (2sfpB) (12) 

and 

| o °° z lf\z) dz = 2A^K 2 {2^B) , (13) 

where K\ is the Mac-Donald's function (modified Bessel's function of the sec- 
ond kind) of order AgI. Based on these equations, one can exactly determine 
the normalization constants A and B. For the case a — 4, the numerically 
determined values of the normalization constants A and B are displayed in 
Table 1. For (3 > 2, we could find the approximate relations 

A w exp (a (3 - 0.1490 a 2 + 1.3689 a + 0.227l) (14) 

and 

B w a (3 + 0.4593 a + 0.9481 . (15) 



4 Data analysis 

We have separately analyzed eight small density intervals in the free low- 
density regime < 20 veh./km and eight density intervals in the congested 
traffic regime > 40 veh./km. After determination of the respective values of (3 
and a from the single vehicle data, we have obtained the fit parameter a by a 
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least square method, i.e. minimization of the error function \ 2 measuring the 
deviation between the theoretical and empirical clearance distributions. The 
predicted Gaussian distributions reproduce the empirical velocity distributions 
very well (see Fig. 1), which is also supported by other empirical and numerical 
studies [6]. The best fit of the netto distance distributions is obtained for the 
integer parameter a — 1 in congested traffic, which is, for example, compatible 
with the intelligent driver model (IDM) and perception-based models [15]. 
Throughout the free traffic regime, we find a good agreement with a — 4, 
corresponding to weak interactions [16] (see Figs. 2, 3). This is not only a 
strong support of studies questioning a uniform behavior of drivers in all 
traffic regimes [3,8]. It also offers an interpretation of the mysterious fractional 
distance-scaling exponent a + 1 ~ 2.8 in classical follow-the-leader models 
[14], which interpolate between the driver behavior in the free and congested 
regimes. 



5 Summary and Discussion 

We have found that it is successful to generalize thermal-equilibrium proper- 
ties of a short-ranged power-law Dyson's gas exposed to a heat reservoir with 
the scaled generalized inverse temperature /3 to steady-state vehicle traffic, 
where j3 is an increasing function of the traffic density p. In the regime of 
congested traffic, this dependency is simply linear: (5{p) = 0.0261p — 0.4785. 
The presented results show that the shape of the interaction potential of ve- 
hicles can be approximated by formula (2) with a = 4 for free traffic and 
a — 1 for congested traffic. The theoretical predictions are further supported 
by the Gaussian distribution of the vehicle velocities in all investigated density 
regimes. 

As it is a hard task to derive analytical relations for the clearance distribu- 
tion (including its normalization constants A and B), we could demonstrate 
the determination of vehicle interaction potentials only for simple functional 
relations. Future advances with this method will hopefully allow to determine 
velocity-dependent interaction forces or functions with turning points, which 
are desireable to reproduce the dynamical behavior in the unstable traffic 
regime realistically as well [18]. Determining interaction potentials in free- 
way traffic does not only contribute to the challenging problem of how to fit 
time headway or distance distributions of vehicles [8,19,2]. It also advances 
the quantitative understanding of human behavior. Normally, it is difficult to 
identify and measure the relevant behavioral variables, and realistic models 
contain a large number of parameters. Here, progress has been made by pow- 
erful methods from statistical physics. Nevertheless, the resulting interaction 
model is not just a physical model. The new picture of two different regimes 
with a = 1 and a = 4 points to adaptive driver behavior to congested and 
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free traffic conditions. The specification of the interaction potential is essential 
for realistic traffic simulations, which are required for the design of efficient 
and reliable traffic optimization measures such as intelligent on-ramp controls, 
driver assistance systems, lane-changing assistants, etc. 
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Table 1 



Numerically determined values of the normalization constants A and B for a = 4. 
The values depend on /?, and therefore, also on the density p (see Table 2). 
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Table 2 



Empirical values of the scaled velocity variation a and the scaled generalized inverse 
temperature (5 for 16 density intervals obtained from single- vehicle data of the Dutch 
two-lane freeway A9, neglecting clearances in front of long vehicles with Zj > 7 m. 
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Scaled Velocity u 

Figure 1. Scaled velocity distributions for eight density regimes in free traffic (above) 
and eight density intervals in congested traffic (below) . The bar diagrams correspond 
to the scaled empirical velocity distributions, while the solid curves correspond to 
the theoretically predicted Gaussian distributions. Note that the mean speeds are 
always scaled to one, while the variances are given by the scaled empirical values 
a 2 (see Table 2). 
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Figure 2. Sum of squared deviations between the empirical and the theoretical netto 
distance distributions for various fit parameters a. The best fit is a = 4 throughout 
the free traffic regime (left) and a = 1 throughout the congested regime (right). 
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Figure 3. Distributions of scaled netto distances (clearances) z = ps for eight den- 
sity intervals in free traffic (above) and eight density regimes in congested traffic 
(below). The bar diagrams represent the scaled empirical distributions, while the 
solid curves correspond to the normalized theoretical distributions with the empiri- 
cally determined values j3 listed in Table 2 and parameter a = 4 for free traffic, but 
a = 1 for congested traffic. 
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